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ABSTRACT 

Photometry and long-slit spectroscopy are presented for a sample of 6 galaxies with a low 



' surface brightness stellar disc and a bulge The characterising parameters of the bulge and disc 

components were derived by means of a two-dimensional photometric decomposition of the 
images of the sample galaxies. Their surface-brightness distribution was assumed to be the 
sum of the contribution of a Sersic bulge and an exponential disc, with each component being 
described by elliptical and concentric isophotes of constant ellipticity and position angle. The 
stellar and ionised-gas kinematics were measured along the major and minor axis in half of the 
^2 ! sample galaxies, whereas the other half was observed only along two diagonal axes. Spectra 

^vq ■ along two diagonal axes were obtained also for one of the objects with major and minor axis 

spectra. The kinematic measurements extend in the disc region out to a surface-brightness 
level fiR ss 24 mag arcsec~ 2 , reaching in all cases the fiat part of the rotation curve. The stellar 
kinematics turns out to be more regular and symmetric than the ionised-gas kinematics, which 
| often shows the presence of non-circular, off-plane, and non-ordered motions. This raises the 

question about the reliability of the use of the ionised gas as the tracer of the circular velocity 
in the modeling of the mass distribution, in particular in the central regions of low surface 
\ brightness galaxies. 

^ , Key words: galaxies: spiral — galaxies: kinematics and dynamics — galaxies: photometry 

— galaxies: structure 
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1 INTRODUCTION 

Low surface brightness (LSB) galaxies are operatively defined as 
galaxies with a central face-on surface brightness fainter than 22.6 
B— magarcsec -2 . Therefore they are more difficult to identify and 
study than their high surface brightne ss (HSB) counterparts. Most 
but not all of them are dwarf galaxies {Schom bert & Bothun|[l988l ; 
ISchombert et alJI 19921 : Umpev et al.lll996t) . Giant LSB galaxies are 
fairly well de scribed as exponential discs with no significant bulge 
components teqmanishin et al.lll983l : iMcGaugh & Bothunlll994 
ISpravberrv et alj 1995t) . They fall within the range of luminosi- 
ties defined by HSB discs. At fixed luminosity, the LSB discs have 



* Based on observations collected at the European Southern Observatory, 
Chile (ESO Prog. 067.B-0283, 069.B-0573, and 070.B-0171). 
f Tables |7l [8l |A6| and lA7l are available only in electronic format 
| E-mail: alessandro.pizzella@unipd.it 



lower central surface brightnesses and larger scale lengths than the 
discs of HSB galaxies, although none of the m has structural p aram- 
eters as extreme as the prototype Malin 1 (Bothu net al.ll987h . It is 
apparent that LSB galaxies exist in a wide range of morphologi- 
cal types. Most LSB galaxies discovered so far are bulgeless, but 
there are als o galaxies with LSB discs that have a significant bulge 
component (Beijersbergen e t al.ll 19991) . 

The amount and distribution of the dark matter (DM) in galax- 
ies are usually determined using the gas rotation curve and through 
the study of an excess of galactic rotation compared to what can be 
induced by the luminous matter, the contribution of which is gener- 
ally estimated by adjusting the mass-to-light ratio until the central 
velocity gradient can be matched. The structure of the DM halo 
is more directly revealed in galaxies where the luminous compo- 
nent gives a nearly zero contribution to the mass budget, because 
this makes mass modeling easier and the derived DM distribution 
less uncertain. This is the case of the giant LSB galaxies. They 
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Table 1. Parameters of the sample galaxies. The columns show the following: (2) morphological classification from RC3, except for ESO-LV 5340200 
(NASA/IPAC Extragalactic Database); (3) heliocentric systemic velocity of the galaxy derived at centre of symmetry of the rotation curve of the gas (this 
paper). The typical error on systemic velocity is AVq = 10 km s _1 ; (4) distance obtained as Vq/ Hq with Hp = 75 km s _1 Mpc -1 and Vo the 
systemic velocity derived from Vq corrected to the CMB reference frame assuming the dipole direction by Fixsen et al. (1996); (5) major-axis position angle 
from ESO-LV; (6) inclinati on derived as cos 2 i = (q 2 — q 2 )/(l — q 2 ). The observed axial ratio q was taken from ESO-LV, except f or ESO-LV 44 60170 
IPalunas & Williamdl200(il) and ESO-LV 2060140 jlVIcGaugh etai]|200lh . The intrinsic flattening q = 0.11 was assumed following iGuthrid ll992h ; (7) 
radius of the 25 B— mag arcsec -2 isophote derived as D25/2 with D25 from ESO-LV; (8) total observed red magnitude from ESO-LV catalog; (9) central 
velocity dispersion of the stellar component (this paper); (10) galaxy circular velocity (this paper); (11) radial extension of the ionized-gas rotation curve in 



units of i?25 (this paper); (12) radial extension of the stellar rotation 


curve in 


units of i?25 (this paper). 
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ESO-LV 1860550 


Sab(r)? 
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40 
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91.7 ±2.0 


235 ± 11 


1.0 


0.6 
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102.7 ±8.2 
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46 
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153.9 ± 7.1 


297 ± 11 


1.4 


0.5 



are DM dominated at all radii when a stellar mass-to-light ratio 
consistent with populatio n synthesis modeli ng is adopted to fit the 
observed rotation curves dSwaters et al . 2000). These results are in 
agreement with the findings of Zwaan et al.l dl995l) . who realized 
that LSB galaxies follow the same Tully-Fisher relation as HSB 
galaxies. They argued that this implies that the mass-to-light ratio 
of LSB galaxies is typically a factor of 2 larger than that of nor- 
mal galaxies of the same total luminosity and morphological type. 
For these reasons, LSB galaxies were considered ideal targets for 
studying the properties of DM halos and for testing if they have 
a central cuspy power-law mass densit y distribution as predicted 
by cold dark matter (CD M) cosmology dNavarro et al.ll 199611 19971 : 
lMooreetal.lll998lll999h . 

Different methods are used to measure the inner slope of the 
DM density profile in LSB galaxies: fitting a mass model which 
takes into ac count the contribution of stars, gas and DM to the ro- 
tation curve (Ivan den Bosch et alj|2000l ; Ivan den Bosch & Swatersl 
l200lMSwaters et al.l2003 fitting a power law to the rotation curve 
dSimon et alj|2003l) . and fitting a power law to the density distribu- 
tion derived by inverting the rotation curve ddeBloketafll200ll ; 
Ide Blok & Bosmall2002l Iswaters et alj|2003ah . For most galaxies 
the inner density slope is unconstrained. Although DM halos with 
constant density cores generally provide better fits, cuspy haloes 
can sometimes be consistent with the data. Systematic errors of 
long-slit (e.g., not-perfect positioning of the slit on the galaxy nu- 
cleus and errors on the position angle of the disc) and radio (e.g., 
beam smearing) observations may lead to an underestimate of the 
inner slope and contribute to the tota l uncertainties on the inter - 
pretation of the observed kinematics ( Swaters et al. 2003a. [20041) 
althou gh only constant density cores can reproduce some observa- 
tions ddeBloketall200ll) . 

High-resolution two-dimensional (2D) velocity fie l ds ob- 
tained with integra l -field spectroscopy 1 Bolatto et alj 20021; 
Swaters et ail l2003bl : iKuzio de Narav et al.l 120061 : IPizzella et ail 



supported (e g iFillmore et alJl986l;lKormendv & Westpfahll 19891 : 



2008) remove most of systematic errors. But, they are not conclu 



sive. The analysis of these data in the inner regions of most of the 
observed galaxies reveals that the gas is not moving onto circu- 
lar orb i ts in a disc coplanar to the stellar one (e.g. , ISwaters et al.l 
l2003bl ; iRhee et al.l |2004| ; lHavashi & Navarro! [2006). It is a com- 
mon phenomenon in the centres of HSB galaxies to o. In fact, 
non-circular (e.g., iGerhardet all I 19891: 1b erma r] |200ll) . pressure- 



iBertola et alTl 19951; ICmzano et alj]l999l). and off-plane (eX 
ICorsini et afll 19991 ; IPignatelli et~ai]|200ll ; ICorsini et al]|2003h gas 
motions are often observed in unbarred HSB galaxies for which 
dynamical modeling allows a direct comparison between the gas 
and circular velocity. Recently. IChristlein & Zarit skv (2008) found 
kinematic anomalies in the outer disk region of nearby disc galax- 
ies. 

This poses a question about the reliability of central mass dis- 
tribution derived in both LSB and HSB galaxies from the analy- 
sis of the gas rotation curves. To address this issue in LSB galax- 
ies a crucial piece of information, which has been missed so far, 
is needed: the stellar kinematics. This is the first of series of pa- 
per aimed at investigating the mass distribution of galaxies with 
LSB discs by means of dynamical models based on both stellar and 
gaseous kinematics. The paper is organized as follows. The galaxy 
sample is presented in §[2] The photometric and spectroscopic ob- 
servations are described in § [3] The structural parameters of the 
bulge and disc of the sample galaxies are derived by analysing their 
2D surface brightness distribution in § [4] The stellar and ionised- 
gas kinematics are measured from the long-slit spectra out to 0.5- 
1.5 optical radii with high signal-to-noise ratio in §|5] The results 
and conclusions are given in §[6] Dynamical modeling of the data 
is deferred to a companion paper (Magorrian et al, in preparation) 



2 GALAXY SAMPLE 

We set out to find LSB discs in spiral galaxies without necessar- 
ily excluding the presence of a bulge. We included both barred and 
unbarred galaxies An initial sample of 10 target galaxies was ini- 
tially selected on the basis of the central disc surface-brightness es- 
timate derived fro m The Surface Photometry C atalog of the ESO- 
Uppsala Galaxies dLauberts & Valentiin|[l989l. ESO-LV). An ex- 
ponential law was fitted to the surface-brightness radial profile at 
large radii, where the light distribution of the galaxy is expected 
to be dominated by the disc contribution. The fit was done using a 
non-linear least-squares minimization me thod. It was based on the 
robust Levenberg-Marquardt method (e.g. iPressetal .1 19921) imple- 
mented by Mor e et al.l dl980l) . The actual computation was done 
using the MPFIT algorithm implemented by C. B. Markwardt un- 
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Table 2. Log of spectroscopic observations of the sample galaxies. The 
columns show the following: (2) observing run; (3) slit position: MJ = ma- 
jor axis, MN = minor axis; (4) slit position angle; (5) number and exposure 
time of the single exposures; (6) total exposure time. 
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X 


2750 


2.3 




2 


MN-30° 


165 


2 


X 


2 750 


1.5 


ESO-LV 2060140 




MJ 


175 


3 


X 


2620 + 1 X 2600 


2.9 




3 


MN 


85 


1 


X 


2600 


0.7 
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Table 4. Log of spectroscopic observations of the kinematical template 
stars. The columns show the following: (2) spectral type from SIMBAD; 
(3) blue apparent magnitude from SIMBAD; (4) observing run; (5) expo- 
sure time; (6) dispersion of the auto-correlation function. 
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Table 3. Log of photometric observations of the sample galaxies. The 
columns show the following: (2) observing run; (3) filter; (4) number and 
exposure time of the single exposures; (5) seeing FWHM of the combined 
image. 
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ESO-LV 5340200 


2 


Gunn z 
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1.1 



der the IDL environment All the sample galaxies are classified as 
spiral galaxies with a morp hological type ranging from Sa to Sm 
(de Vaucou leurs et alJl99ll RC3). The galactic latitudes were cho- 
sen so as to minimize the foreground extinction (|6| > 15°, RC3). 
The inclinations were selected to be larger than 30° (ESO-LV) to 
allow a reliable measurement of the stellar and ionised-gas kine- 
matics. Finally, the galaxies were chosen to be morphologically 
undisturbed and have a clean stellar foreground to not complicate 
the photometric and spectroscopic analysis. 

This sample was subsequently culled to 6 galaxies following a 
scrupulous 2D photometric decomposition of our own images (see 
§ [4j, which led to find that 4 targets were in fact hosting HSB 
discs. The basic properties, photometry, and kinematics of these 
HSB spiral galaxies are given in App. [A] In the rest of the paper 
we discuss only the 6 galaxies w hich host a LSB disc . Our selec- 
tion criteria are similar to those of Beijersbergen et al. ( 1999). The 
result of this selection is that the galaxies in our sample host HSB 
bulges with a central face-on surface-brightness brighter than 22.6 
B— mag arcsec -2 and may have a significant bulge component as 
compared to typical LSB galaxy samples. The presence of a bulge 
is unavoidable when we sample both early and late type spirals . 
The sample galaxies and their basic properties are listed in TableQ] 



1 The updated version of this code is available on 
http://cow.physics. wise. edu~craigm/idl/idl. html 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 Broad-band imaging and long-slit spectroscopy 

The photometric and spectroscopic data of the sample galaxies 
were obtained at the Very Large Telescope (VLT) of the European 
Southern Observatory (ESO) at Paranal Observatory during three 
different runs between 2001 and 2003. The observations were car- 
ried out in service mode on May 05 - July 03, 2001 (run 1), April 
13 - June 03, 2002 (run 2), and October 30, 2002 - February 04, 
2003 (run 3). 

The Focal Reducer Low Dispersion Spectrograph 2 (FORS2) 
mounted the volume-phased holographic grism GRIS.1400V18 
with 1400 grooves mm -1 and the 1 arcsec x 6.8 arcmin slit. In 
run 1 the detector was a Site SI-424A CCD with 2048 x 2048 pix- 
els of 24 x 24 fim 2 . The wavelength range from 4750 to 5 800 A was 
covered with a reciprocal dispersion of 0.5lA pixel -1 and a spatial 
scale of 0.20 arcsec pixel -1 . In run 2 and 3 the detector was a mo- 
saic of two MIT/LL CCDID-20 2k x 4k CCDs with a gap of about 

4 arcsec along the spatial direction. . Each CCD had 2048 x 4096 
pixels of 15 x 15 fim 2 . The wavelength range from 4630 to 5930 
A was covered with a reciprocal dispersion of 0.65A pixel -1 and a 
spatial scale of 0.250 arcsec pixel -1 after a 2 x 2 pixel binning. 

Major and minor-axis spectra were obtained for half of the 
sample galaxies, whereas the other half was observed only along 
two diagonal axes. Spectra along two diagonal axes were obtained 
also for ESO-LV 5340200, which has major and minor axis spectra 
as well. The integration time of the galaxy spectra ranged between 
2400 and 3050 s. The log of the spectroscopic observations as well 
as total integration times and slit position angle are given in Table 2. 
At the beginning of each exposure, the galaxy was centred on the 
slit. The acquisition images were obtained with no filter in run 1. 
The Site CCD efficiency curve roughly corresponded to a broad V 
band. In run 2 and 3 the acquisition images were obtained with the 
Gunn z-band filter (A cff = 9100 A, FWHM = 1305 A). The 
total integration times are listed in Table[3] In run 2 and 3 an offset 
of 10 arcsec along the direction of the slit was applied between 
repeated exposures in order to deal with the gap between the two 
CCDs in the mosaic spectrum. 

In run 1 and 3 a number of spectra of giant stars with spectral 
type ranging from G to K were obtained to be used as templates 
in measuring stellar kinematics. Integration times and spectral type 
the template stars as well as the log of the observations are given 
in Table [4] Additionally, at least one flux standard star per night 
was observed to calibrate the flux of the spectra. The standard cal- 
ibration frames (i.e., biases, darks, and flat-field spectra) as well 
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Figure 1. Surface-brightness radial profile and photometric decomposition 
of the sample galaxies. The radial profiles were obtained by fitting ellipses 
to the isophotes of the image of the galaxy (thick continuous line), model 
bulge (dashed line), model disc (dotted line), and bulge-disc model (dot- 
dashed line). All the profiles are given as a function of the semimajor axes 
of the fitting ellipses. The bulge and disc contributions are shown before 
the convolution with PSF. The surface brightness of the bulge-disc model 
takes into account for PSF. The two thin continuous lines show the ±lcr 
confidence region from the ellipse fit of the galaxy image. 



as the spectra of the comparison arc lamp were taken in the after- 
noon before each observing night. The mean value of the seeing 
FWHM during the spectroscopic exposures as measured by fitting 
a 2D Gaussian to the guide star ranged from 0.8 to 2.7 arcsec. 




[arcsec] 



Figure 2. Comparison between the elliptically-averaged radial profiles of 
surface brightness of ESO-LV 2060140 measured in this work (continuous 
line, with the ±ltr confid ence region indicated by the dotted lines) and by 
Beiiersbergen et "aT Ul999l circles). 



field stars as reference. After checking that their point spread func- 
tions (PSFs) were comparable, the frames were combined to obtain 
a single image. The cosmic rays were identified and removed using 
a sigma-clipping rejection algorithm. 2D Gaussian fits to the field 
stars in the resulting images yielded the final FWHM measurement 
of seeing PSF listed in Table[5] Finally, the mean residual sky level 
was measured in each combined and sky-subtracted image, and the 
error on the sky subtraction was estimated. The median value of the 
residual sky level was determined in a large number of 5 x 5 pixel 
areas. These areas were selected in empty regions of the frames, 
which were free of objects and far from the galaxy to avoid the 
contamination of the light of field stars and galaxies as well as of 
the target galaxy itself. The mean of these median values was zero, 
as expected. For the error in the sky determination we adopted half 
of the difference between the maximum and minimum of the me- 
dian values obtained for the sampled areas. 

For each galaxy, we derived a 'luminosity growth curve' by 
measuring the integrated magnitudes within circular apertures of 
increasing radius by means of the IPvAr0 task ELLIPSE within the 
STSDAS package. Photometric calibration was performed by fit- 
ting this growth curve to that given in R— band for the same circular 
apertures by ESO-LV. The surface-brightness radial profiles of the 
sample galaxies a re shown in Fig.[T] Accord ing to the color gradi- 
ents measured by iBeiiersbergen et al. I d 1999b in a sample of bulge- 
dominated LSB galaxies, we estimated that the accuracy of the ab- 
solute calibration was better t han 0.3 mag arcsec" 2 . The R— band 
radial of ESO-LV 2060140 bv lBeiiersbergen et alj ( fl999h was used 
to test the accuracy of our calibration (Fig. [2}. 



3.2 Reduction of the photometric data 

All the images were reduced using standard MIDAS3 routines. 
First, the images were bias subtracted and then they were corrected 
for pixel-to-pixel intensity variations by using a mean flat field for 
each night. The sky-background level was removed by fitting a 
second-order polynomial to the regions free of sources in the im- 
ages. The different frames of each galaxy were rotated, shifted, and 
aligned to an accuracy of a few hundredths of a pixel using common 

2 MIDAS is developed and maintained by the European Southern Obser- 
vatory. 



3.3 Reduction of the spectroscopic data 

All the spectra were bias subtracted, flat-field corrected, cleaned of 
cosmic rays, corrected for bad columns, and wavelength calibrated 
using standard MIDAS routines. 

The flat-field correction was performed by means of a quartz 
lamp, which were normalized and divided into all the spectra, to 



3 IRAF is distributed by the National Optical Astronomy Observatories 
which are operated by the Association of Universities for Research in As- 
tronomy under cooperative agreement with the National Science Founda- 
tion. 
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Table 5. Photometric parameters of the bulge and disc in the sample galaxies. The columns show the following: (2) effective surface brightness of the bulge; 
(3) effective radius of the bulge; (4) shape parameter of the bulge; (5) axial ratio of the bulge isophotes; (6) position angle of the bulge major axis; (7) central 
surface brightness of the disc; (8) scale length of the disc; (9) axial ratio of the disc isophotes; (10) position angle of the disc major axis; (11) bulge-to-total 
luminosity ratio. 
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ESO-LV 1860550 
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16 ± 0.05 


0.51 ± 0.01 
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111 ± 3 


0.68 


ESO-LV 2060140 


22.27 ± 0.04 


3.1 ± 0.1 
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ESO-LV 2340130 


20.62 ± 0.05 


7.5 ± 0.3 
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72 ± 2 


21.00 ± 0.04 
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ESO-LV 4000370 
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ESO-LV 4880490 
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36.6 ± 0.8 





56 ± 0.01 


88 ± 2 


0.21 


ESO-LV 5340200 


20.44 ± 0.06 


3.0 ± 0.6 
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0.62 ± 0.01 


10 ± 1 
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15.2 ±0.7 





50 ± 0.02 


178 ± 6 


0.56 



correct for pixel-to-pixel sensitivity variations and large-scale il- 
lumination patterns due to slit vignetting. The cosmic rays were 
identified by comparing the photon counts in each pixel with the 
local mean and standard deviation and eliminated by interpolating 
over a suitable value. The residual cosmic rays were eliminated by 
manually editing the spectra. 

The wavelength calibration was performed by means of the 
MIDAS package XLONG. Each spectrum was re-binned using 
the wavelength solution obtained from the corresponding arc-lamp 
spectrum. We checked that the wavelength re-binning had been 
done properly by measuring the difference between the measured 
and predicted wavelengths of about 20 of un-blended arc-lamp 
lines which were distributed over the whole spectral range of a 
wavelength-calibrated spectrum. The resulting rms about the dis- 
persion solution was 0.10 A corresponding to an accuracy in the 
wavelength calibration of 6 km s _1 at 5280 A. The instrumental 
resolution was derived as the mean of the Gaussian FWHMs mea- 
sured for the same un-blended arc-lamp lines we measured for as- 
sessing the accuracy of wavelength calibration. The mean FWHM 
of the arc-lamp lines was 2.04 ± 0.08 A, 2.11 ± 0.07 A, and 
2.10 ± 0.07 A for run 1, 2, and 3, respectively. In run 1 the 1.0- 
arcsec wide slit was projected onto 6 pixels. In run 2 and 3 it was 
projected onto 4 pixels. Therefore, the sampling of the arc-lamp 
lines was lower than in run 1. Nevertheless, FWHM values of the 
different runs are in agreement within the errors. They correspond 
to an instrumental velocity dispersion o"i ns t = 50 km s _1 at 5280 

A. 

All the spectra were correcte d for the misa l ignme nt between 
the CCD mosaic and slit following Bender et al. ( 1994). The spec- 
tra obtained for the same galaxy along the same axis were coadded 
using the centre of the stellar continuum as a reference. This al- 
lowed to improve the signal-to-noise ratio (S/N) of the resulting 
2D spectrum. The accuracy in the wavelength calibration and in- 
strumental velocity dispersion as a function of position along the 
slit were derived from the brightest night-sky emission line in the 
observed spectral range of the coadded 2D spectra of the galaxies. 
This was the [O I] A5577.3 emission line. Its velocity curve and 
velocity dispersi on profile we r e mea sured along the full slit exten- 
sion as done in ICorsini et al. ( 1999). The central wavelength and 
FWHM of the night-sky line were evaluated fitting at each radius a 
Gaussian to the emission line and a straight line to its adjacent con- 
tinuum. Then the values were converted to velocity and velocity 
dispersion, respectively. The velocity curve was fitted by quadratic 
polynomial by assuming the rms of the fit to be the lcr velocity er- 
ror. No velocity gradient was found and the fit rms was 5 km s - . 
The instrumental velocity dispersion derived from the FWHM of 
the [O I] A5577.3 line increased from 46 km s _1 at the slit centre 
to 53 km s _1 at the slit edges. Both the error on the wavelength 
calibration and instrumental velocity dispersion measured from the 



[O I] A5577.3 line in the coadded spectra of the galaxies are con- 
sistent with those derived from the lamp lines in the arc spectra. 

A one-dimensional spectrum was obtained for each kinematic 
template star as well as for each flux standard star. The spectra 
of the kinematic templates were de-redshifted to laboratory wave- 
lengths. In the galaxy and stellar spectra the contribution of the 
sky was determined by interpolating along the outermost 60 arcsec 
at the two edges of the slit, where the target light was negligible, 
and then subtracted. The galaxy spectra always extended less than 
90 arcsec from the centre. A sky subtraction better than 1% was 
achieved. 



4 STRUCTURAL PARAMETERS 

The structural parameters of the sample galaxies were derived by 
applying a 2D bulge-disc photometric deco mposition to their im- 
ages. To this aim the GASP2D algorithm bv lMendez-Abreu et al.l 
d2008l) was used. 

The galaxy surface-brightness distribution was assumed to be 
the sum of the contributions of a bulge and a disc component. The 
Sersic law dSersic 1968) was adopted to describe the surface bright- 
ness of the bulge component. It is parametrized by r e , I E and n 
which are respectively the effective radius, surface brightness at r e , 
and a shape parameter describing the curvature of the radial profile. 
The bulge isophotes are ellipses centred on galaxy centre (xo, yo), 
with constant position angle PAb and constant axial ratio <jb- The 
surface brightness distribution of the disc comp onent was assumed 
to follow an exponential law dFreemanl l 19701) with h and I the 
scale length and central surface brightness, respectively. The disc 
isophotes are ellipses centred on (xo,yo), with constant position 
angle PAd and constant axial ratio g<j . 

Since the fitting algorithm of GASP2D is based on a % 2 min- 
imization, it was important to adopt initial trials for free param- 
eters as close as possible to their actual values. To this aim the 
ellipse-averaged radial profiles of surface brightness, ellipticity, 
an d position angle were analy sed by following the prescriptions 
bv lMendez-Abreu et al.l J2008h . Isophote-fitting with ellipses, after 
masking foreground stars, background galaxies, and residual bad 
columns, was carried out using ELLIPSE. In all cases, first ellipses 
were fitted by allowing their centres to vary. Within the errors, no 
variation in the ellipse centres was found for all the galaxies stud- 
ied in this paper. The final ellipse fits were done at fixed ellipse 
centres. The ellipse-averaged profiles of surface brightness, ellip- 
ticity, and position angle were measured. For all the sample galax- 
ies the ellipse-averaged profiles of surface brightness are shown in 
Fig. Q] The trial values of J c , r c , n, Io, and n were obtained by 
performing a standard photometric decomposition on the ellipse- 
averaged surface-brightness profile. The trial values of PAb and 
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Table 6. Central surface brightness of the LSB discs. The columns show the 
following: (2) effective radius of the galaxy from ESO-LV; (3) mean color 
outside R c from ESO-LV; (4) central blue surface brightness of the disc 
corrected for inclination. 



Galaxy 


R e 


B - R 






(arcsec) 


(mag) 


(mag/arcsec 2 ) 


(1) 


(2) 


(3) 


(4) 


ESO-LV 1860550 


15 


1.23 


22.81 


ESO-LV 2060140 


26 


1.14 


23.08 


ESO-LV 2340130 


18 


1.10 


23.22 


ESO-LV 4000370 


28 


0.97 


23.08 


ESO-LV 4880490 


29 


0.53 


24.88 


ESO-LV 5340200 


11 


0.88 


24.22 



gb were found by interpolating at r c the ellipse-averaged position- 
angle and ellipticity profiles, respectively. The trial values of PAd 
and gd were obtained by fitting a constant to the outer portion of 
the ellipse-averaged position-angle and ellipticity profiles, respec- 
tively. The coordinates of the maximum galaxy surface brightness 
were taken as trial values of the coordinates (xq, yo) of the galaxy 
centre. 

Starting from these initial, trial parameters the bulge-disc 
model of the surface brightness was fitted iteratively by GASP2D 
to the pixels of the galaxy image to derive the photometric param- 
eters of the bulge (J c , r c , n, PAb, and gb) and disc (Iq, h, PAd, and 
gd) and the position of the galaxy centre (xo, yo)- The MPFIT algo- 
rithm based on a robust Levenberg-Marquardt method of non-linear 
least-squares minimization was used in the IDL environment. Each 
image pixel has been weighted according to the variance of its total 
observed photon counts due to the contribution of both galaxy and 
sky, and determined assuming photon noise limitation and taking 
the detector read-out noise into account. The seeing effects were 
taken into account by convolving the model image with a circu- 
lar Gaussian PSF with the FWHM measured using the stars in the 
galaxy image (Table[3]l- The convolution was performed as a prod- 
uct in Fourier domain before the least-squares minimization. 

The parameters derived for the structural components of the 
sample galaxies are collected in Table [5] No correction for galaxy 
inclination was applied. The result of the photometric decomposi- 
tion of the surface brightness distribution of the sample galaxies is 
shown in Fig.[TJ 

The formal errors obtained from the \ 2 minimization method 
are not representative of the real errors in the structural parameters 
dMendez-Abreu etal . 2008). Therefore, the errors given in Table [5] 
were obtained through a series of Monte Carlo simulations. A set of 
3000 images of galaxies with a Sersic bulge and an exponential disc 
was generated. The structural parameters of the artificial galaxies 
were randomly chosen among the following ranges 

0.5 sC r c s; 3 kpc 0.5 ^ g b ^ 0.9 0.5 «C n ^ 6 (1) 

for the bulges, and 

1 s; h <: 6 kpc 0.5 g d s= 0.9 (2) 

for the discs. The artificial galaxies also satisfied the following con- 
ditions 

qd^qb 11 ^ Mr ^ 16 mag. (3) 

The simulated galaxies were assumed to be at a distance of 45 Mpc, 
which corresponds to a scale of 222 pc arcsec -1 . The adopted pixel 
scale, CCD gain, and read-out-noise were respectively 0.25 arcsec 



Table 7. Stellar kinematics of the sample galaxies. The complete table is 
published in the electronic issue. 



Name 


PA 






Av 




Act 


h 3 


Ah 3 


h 4 


Ah 4 


ESO-LV 









km/s 




km/s 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1860550 


75 


-10.02 


102.5 


6.0 


54.8 


2.9 


-0.004 


0.064 


-0.075 


0.133 


1860550 


75 


-16.75 


106.7 


7.4 


63.8 


4.3 


-0.033 


0.051 


-0.080 


0.075 


1860550 


75 


-25.68 


129.4 


8.4 


53.6 


9.4 








00 





NOTE: columns are (1) ESO-LV galaxy number; (2) Position angle on the 
sky; (3) distance from the galaxy center; (4-5) velocity and l-cr uncertain; 
(6-7) velocity dispersion and 1-<t uncertain; (8-9) li3 and l-cr uncertain; 
(10-11) h4 and l-cr uncertain. 

pixel -1 , 2.52 e - ADU -1 , and 3.87 e - to mimic the instrumen- 
tal setup of the photometric observations. Finally, a background 
level and photon noise were added to the artificial images to yield 
a signal-to-noise ratio similar to that of the observed ones. The im- 
ages of artificial galaxies were analysed with GASP2D as if they 
were real. The errors on the fitted parameters were estimated by 
comparing the input pi n and measured p ou t values. The artificial 
and observed galaxies were divided in bins of 0.5 magnitudes. The 
relative errors 1 — Pm/Pout on the parameters of the artificial galax- 
ies were assumed to be normally distributed. In each magnitude 
bin the mean and standard deviation of relative errors of artificial 
galaxies were adopted as the systematic and typical errors for the 
observed galaxies. 

The structural parameters o f ESO-LV 2060140 were ob tained 
from multi-band photometry by Beiiersber gen et at] dl999l) . They 
assumed the bulge and disc to have both an exponential surface 
brightness radial profile and a one-dimensional photometric de- 
composition was applied. The 7?— band bulge (/io,b = 20.19 mag 
arcsec" 2 , hb = 1.62 arcsec) and disc parameters (/xo,d = 21.57 
mag arcsec -2 , hb = 17.75 arcsec) are in good agreement with 
those we derived (Table [5), although we adopted for the bulge 
a Sersic profile with n — 1.64 (fio,b = 20.35 mag arcsec -2 , 
hb — 0.52 arcsec). 

In Table [6] we report the face-on B— band central surface 
brightness of the discs. It was derived from (m (Tab[5} corrected 
for galaxy inclination (Tab. [TJ and was computed by adopting the 
mean (B — R) color given by ESO-LV for radii larger than the 
galaxy effective radius in order to minimize the contamination due 
to the bulge light. 

The six sample galaxies hosting a LSB disc were identified 
among the 10 originally selected galaxies, according to the struc- 
tural parameters derived from the 2D photometric decomposition 
(Tab. |6). All the bulges have a Sersic index n < 2, and few 
of them (ESO-LV 2060140, ESO-LV 4000370) have an appar- 
ent flattening similar to that of the disc (Tab. |6j. According to 
iKormendv & Kennicuttl d2004l) . these disc-like features are the pho- 
tometric signatures of pseudobulges. 



5 STELLAR AND IONISED-GAS KINEMATICS 

5.1 Stellar kinematics and central velocity dispersion 

The stellar kinematics was measured from the galaxy absorption 
features present in the wavelength range between 5050 and 5550 
A, and centred on the Mg line triplet (AA 5164, 517 3, 5184 A) by 
applying the Fo urier Correlation Qu otient method JB ended [l990l 
FCQ) as done in lBenderetal]dl994l) . The spectra were re-binned 
along the dispersion direction to a natural logarithmic scale, and 
along the spatial direction to obtain S/N ^ 20 per resolution el- 
ement. In a few spectra the S/N decreases to 10 at the outermost 
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Figure 3. Examples of the central spectra (continuum line) for three sample 
galaxies. They are compared with the template spectra (dotted line) con- 
volved with their corresponding LOSVDs. The galaxy and template spectra 
have been continuum-subtracted and tapered at the ends with a cosine bell 
function. Residuals are plotted in the upper part of each panel and have false 
zero points for viewing convenience. ESO-LV 2060140 (upper panel) was 
observed in run 3 and fitted with SAO 99192. ESO-LV 4000370 (central 
panel) and ESO-LV 4460170 (see App.fA] lower panel) were observed in 
runs 2 and 1, respectively. They were fitted with SAO 137330. 
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Figure 4. Reliability of the /13 and /14 measurements from Monte Carlo 
simulations. Upper panel: The value of A/i4 obtained for h-j = as a 
function of it* i n and S/N. The values obtained for /14 = are shifted 
from those obtained with /14 = — O.f for viewing convenience. Symbol 
size is proportional to the measured value. Open and filled symbols cor- 
respond to IA/14I smaller and larger than 0.1, respectively. The minimum 
S/N ratios to obtain reliable measurements of I14 for the different cr* - m 
are connected with a continuous line. Lower panel: As in the upper panel 
but for I A/13 1 . 



radii. The galaxy continuum was re moved row-by - row b y fitting a 
fourth to sixth order polynomial as in lBenderetal.1 ( ll994h . The star 
SAO 137330 was adopted as kinematic template for runs 1 and 2, 
and SAO 99192 for run 3. These stars were selected to have the 
narrower auto-correlation function (Table |4j and to minimize the 
mismatch with the galaxy spectra (Fig. [3). 

For each galaxy spectrum the line-of-sight velocity distri- 
bution (LOSVD) was derived along the slit and its moments, 
namely the radial velocity velocity dispersion <r*, and val- 
ues of the coefficients /13 and hi were measured. At each radius, 
they have been derived by fitting the LOSVD with a Gaussian 
plus third- and fourth-order Gauss-Hermite polynomials H3 and 
H4, which describe the asymmetric and s ymmetric deviations of 
the LOSVD from a p ure Gaussian profile ( Ivan der Marel & Franxl 
ll993l;lGerhardlll993r) . The errors on the LOSVD moments were 
derived from photon statistics and CCD read -out noise, calibratin g 
them by Monte Carlo simulations as done bv lGerhard et al. I dl998l) . 
A large number of artificial spectra was generated for each mea- 
sured set of LOSVD moments by first convolving the spectrum of 
the stellar template with the measured LOSVD and then adding 
different realizations of photon and read-out noises. The artificial 
spectra were analysed with FCQ as if they were real. The errors 
on the fitted parameters were estimated by comparing the input 
Pin and measured p ou t values. The relative errors (1 — pi n /p ou t) 
were assumed to be normally distributed, with mean and standard 
deviation corresponding to the systematic and typical error on the 
relevant parameter, respectively. These errors do not take into ac- 
count possible systematic effects due to template mismatch or the 
presence of dust and/or faint emission. 

The stellar velocity dispersion measured at large radii for 



Kinematics of low surface brightness disc galaxies 7 




8 A. Pizzella et al. 




CT -m (km/s) 



Table 8. Ionized-gas kinematics of the sample galaxies. The complete table 
is published in the electronic issue. 
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ESO-LV 
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km/s 


km/s 


km/s 


km/s 


km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


1860550 


75 


-0.75 


23.8 


5.5 


73.7 


4.9 


4.8 


1860550 


75 


0.25 


7.4 


5.7 


75.1 


5.5 


5.4 


1860550 


75 


1.25 


-29.4 


6.2 


61.5 


7.0 


6.8 



Figure 5. Reliability of ct* measurements from Monte Carlo simulations. 
The values of Act* is plotted as a function of ct* only for the artificial 
spectra with reliable measurements of /13 and /14. 



NOTE: columns are (1) ESO-LV galaxy number; (2) Position angle on the 
sky; (3) distance from the galaxy center; (4-5) velocity and 1-ct uncertain; 
(6-8) velocity dispersion and 1-ct lower and upper uncertain following the 
scheme: ct + ^ ct+ 



many sample galaxies is comparable the instrumental velocity dis- 
persion (ct* ~ Uinst ~ 50 km s _1 ). The reliability of such mea- 
surements was assessed by means of a second series of Monte 
Carlo simulations. Artificial spectra were generated by convolv- 
ing the spectrum of the stellar template with a LOSVD with ct* = 
20, 30, 80, 100, 120 km s" 1 , h 3 = and h 4 = 0, -0.1, or 
/13 = 0, — 0.1 and hi = and adopting different values of S/N — 
10, 20, 100. The artificial spectra were measured with FCQ. Fi- 
nally, the Values Of A/13 = /l3,m — /l3,out, A/14 = /l4,in — /U.out 

were computed. The correlation of \Ahz\ and \Ahi\ with ct* and 
S/N is shown in Fig. [4] The measured values of /13 and /14 were 
considered unreliable when \Ahz\ > 0.1 and/or IA/14I > 0.1. 
For this range of ct* and S/N, only and ct* were measured and 
it was assumed /13 = hi — 0. For three of our targets, namelly 
leso-lv 20601401 ESO-LV 4000370 and leso-lv 48804901 h 3 abd hi 
parameters were not measured due to the low signal-to-noise ratio. 
Moreover, the stellar velocity dispersion was systematically under- 
estimated for ct* ^ 30 km s _1 . This was found by measuring 
Act* = CT*.i n — CT* iOUt in artificial spectra with IA/13I ^ 0.1 and 
I A/14 1 ^ 0.1 (Fig. [5}. Such a small value of ct* is far below the 
velocity dispersions measured for the sample galaxies indicating 
that these measurements are genuine and are not strong affected by 
instrumental broadening. This is not the case of the small velocity 
dispersions (~ 25 km s" 1 ) measured for the faint foreground stars 
(S/N ~ 20) which happened to fall into the slit when some of 
the galaxy spectra w ere obtained. These results are consistent with 
previous findings bv lJoseph et al. I j200lh : |Pmknev et alj J2003h . 

The measured stellar kinematics of all the sample galaxies are 
reported in Tableland plotted in Fig.[6] The central velocity disper- 
sion of the stellar component ct c was derived by extrapolating the 
velocity dispersion radial profile to r = arcsec. This was done by 
fitting the eight innermost data points with an empirical function. 
The best-fitting function was selected to be either exponential, or 
Gaussian, or co nstant. No a pertur e cor rection was applied to ct c , 
as discussed by iBaes et all 020031) and IPizzella et alj J20041) . The 
derived values for ct c of our sample galaxies are given Table[T] 

5.2 Ionised-gas kinematics, circular velocity, and minor-axis 
rotation 

The ionised-gas kinematics was measured by the simultaneous 
Gaussian fit of the H/3 and [O III]AA4959, 5007 emission lines 
present in the galaxy spectra. The galaxy continuum was removed 
from the spectra as done for measuring the stellar kinematics. A 
Gaussian was fitted to each emission line, assuming they have the 
same line-of-sight velocity v s , and velocity dispersion a g . An addi- 
tional absorption Gaussian was considered in the fitting function to 



take into account the presence of the H/3 absorption line. Far from 
the galaxy centre, adjacent spectral rows were average to increase 
the signal-to-noise ratio (S/N > 10) of both the em ission lines. 
The fitting routine developed by IPizzella et alj J2004l) is based on 
the MPFIT algorithm to perform the non-linear least-squares mini- 
mization in the IDL environment. 

We checked for some galaxies that the error in the kinematic 
parameter determination derived by Monte Carlo simulations did 
not differ significantly from the formal errors given as output by 
the least-squares fitting routine. We therefore decided to assume 
the latter as the errors on the gas kinematics. For each galaxy 
we derived the heliocentric systemic velocity, Vq, as the veloc- 
ity of the centre of symmetry of the ionised-gas rotation curve 
(Table [T). The measured ionised-gas kinematics of all the sam- 
ple galaxies are reported in Table [8] and plotted in Fig. [6] The 
ionised-gas kinematics along the major axis of some of the sam- 
ple galaxies have been measured by other authors. A comparison 
with these data sets was performed to assess the accuracy and re- 
liability of our measurements (Fig. [Til . The Ha ro tation curves ob- 
tained for ES O-LV 2060140 dMcGaugh et alj200ll Fig.|7Ji), ESO- 
LV 2340130 dMathewson etalj|l992t Fig. 17b). ESO-LV 4880490 
(McGaugh et alj200ll Fig.rTb). and ESO-LV 5140100 (see App.[A] 
iMathewson et alj 1992k Fig.[7j)), are in agreement within the errors 
with those measured by us after applying a small offset in their sys- 
temic velocities (AV ^ 30 km s~ ). As a general comment, we 
can notice that our data are characterized by a finer spatial sam- 
pling. This allows to reveal kinematic details that were not evident 
in previous observations. The extension of the kinematic data are 
comparable. It has to be noted that our measurements are based on 
the H/3 and [O III] emission lines, while the data available in the 
literature are based on the Ha emission line. 

For all the sample galaxies, the ionised-gas rotation curve 
measured along the major axis (or along diagonal axis which was 
closer to the major axis) was obtained by folding the observed line- 
of-sight velocities around the galaxy centre and systemic velocity 
after averaging the contiguous data points. The curve was depro- 
jected by taking into account the inclination of the disc (Table [TJ 
and the misalignment between the observed axis and the position 
angle of the galaxy major axis (Table 2). The galaxy circular ve- 
locity V c was derived by averaging the ou termost values of the flat 
portion of deprojected rotation curve (see IPizzella et alj d2005l) for 
a detailed descrption). Their V c are listed in Table Q] For five out 
of six sample galaxies the ionised-gas kinematics were measured 
either along the minor axis (ESO-LV 2060140, ESO-LV 4880490, 
ESO-LV 5340200) or along a diagonal axis close to the minor one 
(by 30° and 16° for ESO-LV 1860550 and ESO-LV 2340130, re- 
spectively). The corresponding rotation curves are either irregu- 
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Figure 6. Kinematic parameters for stars (circles) and ionized-gas (crosses) measured along the observed axes of the sample galaxies. The radial profiles of the 
line-of-sight velocity (after the subtraction of systemic velocity and with no correction for inclination), velocity dispersion (corrected for instrumental velocity 
dispersion), third and fourth order coefficient of the Gauss-Hermite decomposition of the stellar LOSVD are shown (from bottom to top). For each panel the 
galaxy name, location and position angle of the slit are given. 



lar and asymmetric (ESO-LV 1860550, ESO-LV 2060140, ESO- 
LV 4880490) or characterized by a well-defined velocity gradient 
(ESO-LV 2340130, ESO-LV 5340200). Significant velocity values 
are unexpected along the galaxy minor axis if the gas is moving in 
circular orbits in a disc coplanar to the stellar one. The asymmetry 
of rotation curve of ESO-LV 1860550 (PA = 165°) and the rever- 
sal of the velocities measured for ESO-LV 2340130 (PA = 98°) 
can not be attributed to the small misalignment (APA = 30° and 
16° ) between the slit and the galaxy minor axis. They are sug- 
gestive of a more complex gas distribution. As a consequence, the 
ionised-gas velocity measured along the major axis is not tracing 
the circular velocity in the inner regions of the galaxy. To properly 
address this issue it is crucial to map the entire velocity field of gas. 
This is the case of ESO-LV 1 860550, ESO-LV 4000370, and ESO- 
LV 5340200. The 2D velocity field of their gaseou s component 
has been recently measured by IPizzella et al.l j2008t) . They found 
the presence of non-ordered gas motions in ESO-LV 1860550, of a 
kinematically-decoupled component in ESO-LV 4000370, and of a 
misalignement between the photometric and kinematic major axis 
of ESO-LV 5340200. 



6 SUMMARY AND CONCLUSIONS 

We have presented and investigated the structural parameters and 
kinematic properties of 6 galaxies with a LSB disc, for which 
we obtained deep photometric and spectroscopic VLT observa- 
ti ons. Initially chosen follow ing selection criteria similar to those 
of Beii ersbergen et al.l ( 1999) and using published photometric pro- 
files, we have subsequently confirmed the presence of a LSB disc 
in our sample galaxies by carrying out a careful 2D photometric 
bulge-disc decomposition. Our sample galaxies differs as compared 
to typical LSB galaxy samples for the possible prsence of a signif- 
icant bulge component 

Our photometric analysis further revealed that our LSB galax- 
ies have bulges that can are well described by a Sersic profile with 
an index n < 2, with few of them displaying an apparent flat- 
te ning rather similar to th a t of th eir surrounding disc. According 
to lKormendv & Kennicuttl J2004l) . these disc-like features are the 
photometric signature of a pseudobulge. 

We obtained spectroscopic observations along the major and 
minor axis in half of our sample galaxies (ESO-LV 2060140, ESO- 
LV 4880490 and ESO-LV 5340200,), while we observed the other 
half along two diagonal axes. Spectra along two diagonal directions 
were obtained also for ESO-LV 5340200. In all targets, the depth of 
our data allowed us to measure the stellar and ionised-gas kinemat- 
ics out to a surface-brightness level (ir ~ 24 mag arcsec -2 , com- 
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Figure 6. Continued 



fortably reaching the outer disc regions where the rotation curve 
flattens. 

We derived central values for the stellar velocity dispersion 
of (j c and measured the circular velocity in the outskirts of our 
sample gal axy V c by follo wing the flat portion of the gas rota- 
tion curve. iFerraresd J2002 ) found a correlation between a c and V c 
and interpreted it as evidence for a connection between the mass 
of the c entral supermassive black hole and that of the DM halo 
(see also lBaes et alJl2003l ; lBuvle et alj|2006h . Using the a c and V c 



measurements presented in this paper, IPizzella et alj |2005) found 
that galaxies hosting LSB discs follow a different V c — a c relation 
with respect to their HSB counterparts. In fact, they show either 
higher V c or lower a c . This behaviour has been recently explained 
in terms of morphology or, equivalentl y, total light concentration 
bv lCourteau eTal] d2007t) andE3 d2007l) . 

We have shown that the vast majority of our sample galaxies 
show non-zero gas velocities along or close to their minor axes. 
Significant velocities are unexpected along the galaxy minor axis if 
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the gas is moving in circular orbits in a disc coplanar to the stellar 
one. They can be explained as due to non-ordered motions or by 
invoking a warp and/or non-circular motions of the gaseous com- 
ponent. As a consequence, the ionised-gas velocity measured along 
the major axis is not tracing the circular velocity in the inner regions 
of the galaxy. Deviations from perfect circular gas motions are by 
no means limited to the minor axis direction in our sample galaxies. 
Along the major axis of our targets (or the closest diagonal direc- 
tion) the ionised-gas velocity curve is characterized by an higher 
degree of asymmetry and irregularity with respect to the stellar one. 
As a consequence, to derive the gas rotation curve it is often nec- 
essary to smooth the data over large radial bins before folding the 
gas velocities around the centre of symmetry. This is not the case 
for the stars, since their measured velocities fold quite well and 
do not show sharp wiggles and bumps. The rotation curves of the 
ionised gas and stars are plotted in Fig.[8]for all the sample galaxies. 
Data were folded around the centre of symmetry and deprojected 
by taking into account the inclination of the disc (Table[]} and the 
misalignment between the observed axis and the position angle of 



the galaxy major axis (Table 2). Globally the gas rotation curves 
of the sample galaxies are more asymmetric, irregular, and affected 
by scatter than the corresponding stellar ones. The presence of a 
more disturbed gas motions and a regular stellar kinematics is par- 
ticularly clear in ESO-LV 1860550 and ESO-LV 2340130. Never- 
theless, locally both gas and stellar rotation can both appear pretty 
smooth, in particular toward the centre (e.g., ESO-LV 2340130 and 
ESO-LV 5340200). On the other hand, the most obvious scatter of 
the gas data points at larger radii is instrinsic and neither due to 
low S/N nor to the underlying H/3 asborption line, since this is ob- 
served only at small radii (where the kinematics is constrained also 
by [O III] emission doublet). 

While it is easy to spot non-circular, off-plane, and non- 
ordered velocities along the minor axis, their presence along 
major axis may remain undetected. This can affect the correct 
estimate of the central velocity gradient, posing severe limitations 
to the use of the ionised gas as the tracer of the circular velocity 
to derive the mass content and distribution in the central region 
of galaxies. To address this issue, it is crucial to measure the 
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kinematics along different axes and ideally to have integral-field 
data. The 2D velocity field of the ionised gas in LSB galax- 
ies has been mapped only recently by means of integral-field 
spectroscopy jKuzio de Narav et al] 120061 : ICoccato et alj 1200 8t 
IPizzella et aUl2008h . IPizzella et al] j2008) reported the detection 
in ESO-LV 1890070 and ESO-LV 4000370 of strong non-ordered 
motions of the ionised gas, which prevent to disentangle between 
core and cuspy radial profiles of the DM mass density. 

The main conclusion of this work is that the stellar kinemat- 
ics of galaxies with LSB discs is less affected than the ionised-gas 
kinematics one by large-scale asymmetries and small-scale irreg- 
ularities. The gas velocities have been found to deviate from the 
from the circular velocity in particular and more often in the inner 
regions of our sample galaxies. Unless the regularity of the ionised- 
gas motions can be establish with integral-field data, our findings 
strongly suggest that stars are more suitable to trace the mass dis- 
tribution of galaxies with a LSB disc, in particular when the inner 
slope of the dark matter halo is the central question that one wishes 
to answer. 

Building a complete dynamical model to address the question 
of the mass distribution of our sample galaxies is beyond the scope 
of this paper, and will be presented in a companion paper (Magor- 
rian et al, in preparation) 
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APPENDIX A: HSB GALAXIES 

ESO-LV 1890070, ESO-LV 4460170, ESO-LV 4500200, andESO- 
LV 5140100 (Tab. I Alb were originally selected according to the 
criteria given in §|2] Their photometric (Tab, |A3t and spectroscopic 
(Tab, IA2t observations were carried out and data analysed together 
with the sample galaxies, as discussed in the paper. 

In spite of the selection criteria, ESO-LV 1890070, ESO- 
LV 4460170, ESO-LV 4500200, and ESO-LV 5140100 do not 
host a LSB disc, according to the structural parameters derived 
from the 2D photometric decomposition (Tab. lA4l (. The central sur- 
face brightness of the discs derived from ESO-LV photometry was 
overestimated with respect to that from the available CCD images 
(Tab. |A5I >. This is due to the limited radial extension of ESO-LV 
data and to the shallow surface-brightness radial profile observed 
at the largest measured radii. The flattening of the light profile is 
actually due to the presence of a bar in ESO-LV 1890070, ESO- 
LV 4460170, and ESO-LV 4500200 (as found by subtracting the 
photometric model from the observed image of the galaxy), and 
by the downbending of the disc radial profile in ESO-LV 5140100 
(Fig.fAD. 

The stellar and ionised-gas kinematics of the HSB galaxies is 
plotted in Fig. IA2l and given in Tables lA6l and [ATI These data al- 
lowed to derive the bulge velocity dispersio n g mc and disc circu lar 
velocity V c (Tab. lAlt . which were used by Pizz ella et al.l d2005f) to 
study the V c — <J C relationship in HSB and LSB galaxies. The rota- 
tion curves of the ionised gas and stars are plotted in Fig. [A3] 

All the galaxies show non-zero velocities for the ionised gas 
and stars along the minor axes. The minor-axis velocity gradi- 
ents measured for ESO-LV 1890070, ESO-LV 4460170, and ESO- 
LV 4500200 can be explained as due to the non-circular mo- 
tions induced by the bar. Moreover ESO-LV 1890070 and ESO- 
LV 5140100 host a dynamically cold component in their nucleus. 
The a+ radial profile of ESO-LV 1890070 along both the major and 
minor axis are characterize d by a central minimum . A similar phe- 
nomenon was observed by Emsellem et al. (2001) in a number of 
barred galaxies and interpreted as the signature of a nuclear stellar 
disc or bar. The presence of such a fast-rotating nuclear component 
in ESO-LV 1890070 is also supported by the sharp stellar veloc- 
ity gradient (~ 300 km s _1 kpc" 1 ) measured along the galaxy 
major axis. A central drop was observed in the cr* radial profile of 
ESO-LV 5140100 too. But this is an unbarred galaxy according to 
the results of the photometric decomposition. Therefore, the lower 
central value of <r* has to be attributed to the small, low-luminosity, 
and exponential bulge of the galaxy. 

This paper has been typeset from a TgX/ LTr-X file prepared by the 
author. 
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Table Al. As in Tab.Qjbut for the HSB galaxies. The inclination of ESO-LV 4460170 is from lPalunas & Williams! feOOd) . 
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Table A4. As in Tab.[5]but for the HSB galaxies. 



Galaxy fie r e n q h PA b ^0 h 9d PA d 

(mag/arcsec^ ) (arcsec) (°) (mag/arcsee^ ) (arcsec) (°) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 

ESO-LV 1890070 19.24 ± 0.03 3.4 ±0.1 0.99 ± 0.02 0.71 ±0.01 178 ± 4 20.32 ± 0.03 31.1 ± 0.8 0.46 ±0.01 3±1 

ESO-LV 4460170 18.97 ± 0.03 3.0 ±0.1 1.23 ± 0.03 0.85 ± 0.02 177 ± 4 20.11 ±0.03 19.4 ± 0.5 0.51 ±0.01 1±1 

ESO-LV 4500200 18.01 ± 0.03 3.3 ±0.1 1.79 ± 0.04 0.70 ±0.01 179 ± 4 19.03 ± 0.03 17.3 ± 0.4 0.67 ±0.01 4±1 

ESO-LV5140100 20.19 ± 0.03 2.2 ±0.1 0.93 ± 0.02 0.91 ± 0.02 171 ± 4 20.45 ± 0.03 27.1 ± 0.7 0.78 ± 0.02 175 ± 5 



Table A2. As in Tab.[2]but for the HSB galaxies. 



ESO-LV 1890070 



ESO-LV 4460170 



Galaxy Run 


Position 


PA 


Single E\p. Time 


Tot. Exp. Time 






(°) 


(s) 


(h) 


(1) (2) 


(3) 


(4) 


(5) 


(6) 


ESO-LV 1890070 1 


MJ 


18 3 X 


2400 


2.0 




MN 


108 2 X 


2900 


1.6 


ESO-LV 4460170 1 


MJ 


155 3 X 


2400 


2.0 




MN 


65 2 X 


3000 


1.7 


ESO-LV 4500200 1 


MJ 


154 3 X 


2400 


2.0 




MN 


64 3 X 


2400 


2.0 


ESO-LV 5140100 1 


MJ 


50 3 X 


2520 ± 2 X 2400 


4.8 




MN 


140 1 X 


2400 


0.7 


Table A3. As in Tab.[3]but for the HSB galaxies. 


Galaxy 


Run 


Filter Single Exp. Time 


FWHM 








(s) 


(arcsec) 


(i) 


(2) 


(3) 


(4) 


(5) 


ESO-LV 1890070 


1 


none 


4 X 20 


0.8 


ESO-LV 4460170 


1 


none 


4 X 20 


1.4 


ESO-LV 4500200 


1 


none 


2 X 20 


1.1 


ESO-LV 5140100 


1 


none 


4 X 20 


1.1 


Table A5. As in Tab.|6]but for the HSB galaxies. 


Galaxy 


R c 


B — R 








(arcsec) 


(mag) 


(mag/arcsec 2 ) 




(1) 


(2) 


(3) 


(4) 




ESO-LV 1890070 


48 


0.89 


21.67 




ESO-LV 4460170 


24 


1.24 


21.92 




ESO-LV 4500200 


23 


1.68 


20.87 




ESO-LV 5140100 


42 


1.16 


21.84 






50 100 150 200 30 40 60 80 100 
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30 40 60 80 100 20 40 60 80 100 120 
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Figure Al. As in Fig.[TJbut for the HSB galaxies. 



Table A6. Stellar kinematics of the HSB galaxies. The complete table is 
published in the electronic issue. 



Name 


PA 






Av 




Act 


h 3 


Ah :i 


h 4 


Ah 4 


ESO-LV 







km/s 


km/s 


km/s 


km/s 








(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1890070 


108 


-94.50 


-139.4 


25.5 


117.7 


22.8 


0.00 


0.00 


0.00 


0.00 


1890070 


108 


-70.64 


-131.4 


8.5 


63.5 


7.0 


-0.01 


0.09 


-0.07 


0.13 


1890070 


108 


-62.55 


-127.8 


6.4 


58.4 


4.6 


0.00 


0.07 


-1106 


0.12 



NOTE: columns are (1) ESO-LV galaxy number; (2) Position angle on the 
sky; (3) distance from the galaxy center; (4-5) velocity and l-cr uncertain; 
(6-7) velocity dispersion and l-cr uncertain; (8-9) li3 and l-cr uncertain; 
(10-11) h4 and l-cr uncertain. 
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Figure A3. As in Fig.[8]but for the HSB galaxies. 



Table A7. Ionized-gas kinematics of the HSB galaxies. The complete table 
is published in the electronic issue. 



Name 


PA 


r 


V 


Av 


a 


Act- 


A(T+ 


ESO-LV 


O 




km/s 


km/s 


km/s 


km/s 


km/s 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


1890070 


108 


-125.30 


-151.6 


6.2 


13.6 


13.6 


92.2 


1890070 


108 


-123.30 


-135.3 


5.8 


22.7 


7.0 


7.0 


1890070 


108 


-121.30 


-139.2 


5.4 


25.9 


4.5 


4.5 



NOTE: columns are (1) ESO-LV galaxy number; (2) Position angle on the 
sky; (3) distance from the galaxy center; (4-5) velocity and l-cr uncertain; 
(6-8) velocity dispersion and l-cr lower and upper uncertain following the 
scheme: cr + ^ CT+ 



